Search for: All records

Symmetry-adapted perturbation theory (SAPT) is an ab initio approach that directly computes noncovalent interaction energies in terms of electrostatics, exchange repulsion, induction/polarization, and London dispersion components. Due to its high computational scaling, routine applications of even the lowest order of SAPT are typically limited to a few hundred atoms. To address this limitation, we report here the addition of electrostatic embedding to the SAPT (EE-SAPT) and ISAPT (EE-ISAPT) methods. We illustrate the embedding scheme using water trimer as a prototype example. Then, we show that EE-SAPT/EE-ISAPT can be applied for efficiently and accurately computing noncovalent interactions in large systems, including solvated dimers and protein–ligand systems. In the latter application, particular care must be taken to properly handle the quantum mechanics/molecular mechanics boundary when it cuts covalent bonds. We investigate various schemes for handling charges near this boundary and demonstrate which are most effective in the context of charge-embedded SAPT. 

Abstract The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA. 

Abstract Cyanuric acid (CA), a triazine heterocycle, is extensively utilized for noncovalent self‐assembly. The association between poly(adenine) and CA into micron‐length fibers was a remarkable observation made by Sleiman and co‐workers, who proposed that adenine and CA adopt a hexameric rosette configuration in analogy with previously reported structures for CA assemblies. However, recent experimental observations from the Krishnamurthy group led to a reevaluation of the hexameric rosette model, wherein they have proposed a hydrogen‐bonded helicene model as an alternative. Our molecular dynamics simulations show that the hexad model is indeed unlikely and that this novel noncovalent helicene geometry, where the adenine and CA bases form an extended helical hydrogen‐bond network across the system, is a more probable structural motif. The existence of noncovalent helicene compounds may have wide‐ranging applications in DNA nanotechnology and helicene chemistry.